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1. INTRODUCTION 

The problem of the behavior of a real quan- 
tum system is complicated in many cases by the 
fact that it is coupled to the environment de- 
grees of freedom. Therefore, the adequate de- 
scription of the non-equilibrium time-dependent 
properties of quantum systems in contact with a 
thermal reservoir and in external time-dependent 
fields is an interesting and important problem of 
modern theoretical physics. In particular, this 
description is necessary for studying stochastic 
processes in mathematical physics [1 121,31 lj and 
for studying the effective potentials at finite tem- 
peratures in particle physics [5]. In condensed 
matter physics examples of this type of problems 
are two-state systems coupled to a dissipative en- 
vironment 6], spin systems in a heat bath en- 
vironment [217181911011 H12ll"3| . superconductors 
in contact with a heat bath |14ll5IH)ll7ll8ll9l20j 
and many others. 

The standard method for studying quantum 
systems out of equilibrium is the closed-time-path 
Green function formalism |21I22| . For spin sys- 
tems in general, this method can't be applied di- 
rectly, since the commutators of arbitrary spin 
operators S are not c-numbers, and for this case a 
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standard Wick theorem does not exist [2U • How- 
ever, for the particular case of S = ^, when the 
operators anti-commute, the corresponding path 
integral technique has been developed |7I24I25| . 
This technique is based essentially on standard 
field theory methods with the introduction of Ma- 
jorana fermions. In this paper we study the non- 
equilibrium properties of the S — 1/2 ferromag- 
netic Heisenberg model in contact with a phonon 
heat bath by means of this method (the Heisen- 
berg model in equilibrium was studied with this 
technique in .8 26 , for example). 

There are several reasons for studying spin sys- 
tems out of equilibrium. For example, recently 
discovered big magnetic molecules |27|28l29l30| 
show relaxation times which are extremely large 
(see also where the experimental results for 
the case of external magnetic field are described). 
Theoretical studies of this phenomenon were per- 
formed in |9I1()I12| . but more quantitative analy- 
sis is still lacking. 

The behavior of a spin system in the presence 
of an external magnetic field flip H — > — H is 
also an interesting problem. In particular, the 
question whether there exists a critical value of 
the magnetic field below which the system does 
not relax to equilibrium and is in a quasi-periodic 
regime, as it takes place in the case of a 4> -theory 
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|32I33| . is interesting. 

In this paper the questions mentioned above 
are studied in the case of the Heisenberg model 
coupled to a phonon heat bath and in an exter- 
nal time-dependent magnetic field. We consider 
the interacting spin system in a magnetic field on 
the mean field level, and the phonon heat bath 
is studied as a perturbation up to the second or- 
der in the spin-phonon coupling. We introduce 
the formalism specific to the spin problem and 
apply it to several examples of a single spin in a 
magnetic field and then take the leading approx- 
imation to study interacting spins. 

2. FORMALISM 

The Heisenberg Hamiltonian with a spin- 
phonon coupling in an external magnetic field can 
be written as: 

H = — — JijSiSj — HiSi 

ij i 

+ + \ + m^4) (1) 

qi qi H 

where Jjj > is the nearest neighbor ferromag- 
netic coupling, i, j are site coordinates, q is a bath 
mode with coordinate X q i, momentum p q i and 
frequency u> q i and c q i is the spin-phonon coupling. 

The normalized generating functional of the 
system is 



Z 



(2) 

where H s and Hi, are the spin and bath parts of 
the Hamiltonian, respectively. Initially the spins 
and the phonons are decoupled. The spin system 
is at a temperature Tj = 4- and the phonon bath 
is at a temperature Tf = At time ti the spin- 
phonon coupling is turned on and it is assumed 
that the spin system will reach thermodynamic 
equilibrium with the phonons. In the closed-time 




Figure 1. Integration contour for the time vari- 
able. The direction is U — > tf — > ti — > i, ti — 



path Green function formalism the numerator of 
the generating functional can be written in this 
case in terms of the path integral: 
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(3) 



where D is the phonon propagator and the time 
integration contour C is presented in Fig.l. 

In (J3J) a Majorana fermion operator Q is intro- 
duced. In the case of the spin ^ the following 
relation holds [2317] : 



Si — — Ci X Ci 



(4) 
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After the integration over the phonon modes, the 
numerator of the generating functional (J3J) be- 
comes 



D(exp 



i I dt 
'c 



- — ^ Si Jij Sj Hi 



and the trace of the heath bath in the denomina- 
tor is cancelled. The presence of the phonon bath 
environment is given by the Feynman and Vernon 
influence functional |1I2| . Its effect is to change 
the spin-spin coupling to an effective coupling 

Mt - O -» - + iA v(t - O- ( 5 ) 

The function Ay-(i — t') is 

Aij(t-t') = a l] {t-t')6{t-t') + a* :j (t-t')6(t' -t), 

e iu q (t—t') e -iu) q (t-t') 



aij(t-t') = 



C qi C qj 

2m q u)q 



1 1 



-f3 f u 



and it contains all the effective properties of the 
bath environment. It can be approximated in the 
standard way |2J: 



a i3 (t-t') 



S tj A 



a iui(t-t') 



dlOLU S 



— iul(t—t') 



o/u 



1 1 



-0fO> 



where A 7 s and ui c are effective phonon density of 
states parameters . In particular, A is an effective 
spin-phonon coupling, lu c is the phonon energy 
cut-off, and s describes the "subohmic" , "ohmic" 
and "superohmic" cases, when < s < 1, s = 1 
and s > 1, respectively. Below we shall consider 
separately these different cases. 

To decouple the interaction between the spins 
let us make the Hubbard- Stratonovich transfor- 
mation: 



DC,D4> 
det(2irJ) 



exp 



-t I dt 
lc 



*3 



iJij 1 ^! - E fa^i ~ /;HjSi 



and with the shift 



H, 



Z = 



exp [ — i I dt 

det(2wJ) \ Jc 



In the case of the substitution (J2J the effective 
fcrmion Hamiltonian is 



(7) 



Heff = 



& = § E ' 

ilran 



I hnn f-m >n 



where l,m,n are vector components. The ac- 
tion in the exponent of (J7J is quadratic on the 
^■-variables. Therefore it can be formally inte- 
grated, 



Z = 



D(f) 



xe 

where 
K(Si) 



'det(2TtJ) 

lo^^-^^-^YlK^A, ( 8 ) 



D( exp 



dt 



c 



(9) 



is a function of the order parameter <f>. The func- 
tional K((j>i) is the Helmholtz free energy for a 
single spin in an external effective field fa. 

If the external field is uniform in space and 
time, the propagator of the £-fields 

G${t,t') = -i<fQ{t)C]{t')> (10) 

can be obtained analytically. In the case of an 
external magnetic field H = (0,0, H z ), this can 
be done most conventionally in the spherical basis 
with unit vectors a± = ^(s x ± ie y ), do = e z [B]. 
It becomes diagonal and satisfies the equation 

(i± aH *)G a (t-t') = 5(t-t'), 
dt 

where a — ±, denotes the spherical components. 
The times t and t' take values on the contour C, 
and therefore G can be presented as 



(6 ) G a (t-t') = G>(t-t')9 c (t-t') 

+G<(t-t')e c (t'-t), 



(ii) 
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(t-tj). 



Figure 2. Time-dependence of one-loop magneti- 
zation for one spin after magnetic field flipping at 
different values of s and with A = 0.025, H/u> c = 
0.1,T i /w c = 0.1,T / /w c = 0.02. 



where 9c it' ~ t) is the theta- function on the 
contour C. 

The solution of (|f 0|l in the case of a time- 
independent magnetic field H with the Kubo- 
Martin-Schwinger boundary conditions |34) 

G>(t-i^-l/) = -G<{t-l/) 



is 



G>(t-t') = -i- 



^-iuH z (t-t' 



G<(t-t') = i 



1 



1 + gftaff* 



-iaH z (t-t') 



\ (12) 



(13) 



However, in general, since the field may be both 
time and space dependent there is no simple an- 
alytical solution. Therefore, in practice, the £- 
integration in can be performed using the dis- 
crete path integral approximation, with the time 



contour divided on a finite number N of time in- 
tervals and Eq.@ written as: 



Ki4>) = 
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and and 1 are diagonal 3x3 zero- and unit- 
matrices, respectively. The discrete inverse prop- 
agator in the presence of a general magnetic field 
<b is then 



G~ 



with 



-2i 



( c - 1 i . 
i c -i 
-iic.. 



V i 



l l 
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- i 
'c J 



At 



(15) 
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The integration over the ^-fields in @ gives 
1 f D$ 



Z 



23N/2 



det (2tt J) 



x exp 



dt^^i-HiW/^j-H^ 

+~ Tr ln iG_1 (^)J = (17) 

where G(cj>i) is calculated numerically. The loop 
expansion on the field tj> in the thermodynamic 
potential can be performed to take into account 
the fluctuations of the order parameter • 

Minimization of the thermodynamic potential 
with respect to <f> gives the saddle-point equation 
for the order parameter: 

$f\t) =Hi{t)+ f dt'J l0 it,t')Mf\t'), (18) 
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Figure 3. Time-dependence of one-spin magne- 
tization after lowering temperature for the cases 
of magnetic field nipping and constant magnetic 
field at s = 0.5, A = 0.025, H/uj c = 0.1, T,/lu c = 
0.1 and Tf/w c = 0.02. 



where the mean-field magnetization Mj(t') is 



l. 



~e lsr G%(t',t'). 



Taking the time limit properly, one gets 

M? \t) = y i »r(G%(t + ,t) + G° r j (t-,t)). (19) 

The self-consistent system of the equations <|18|) . 
lfT^|l together with l(T5|l - (fTfi|l can be solved nu- 
merically, giving the ^-dependence of the Green 
function (|10(l . which allows to analyze the mag- 
netic properties of the system. In this paper we 
consider the spatially homogeneous magnetic field 
case, so the Green's function for the ^-field can be 



found analytically in some cases of the field time 
dependence. 

The magnetization M^ x \t) calculated to the 
order of one- loop in the boson fluctuations is |26j : 

5(-0FW) 



M* (1) {t) 



where 



SHHt) ' 



(20) 



is the free energy calculated to the order of one 
loop in the effective RPA interaction D: 



D^-\t,t') = J^- 1 S(t,t')-ix P s T(t,t'), 
with the free magnetic susceptibility 

x7(*,f) 

_ _}_ £ lhh e m mi m 2G hp ^ t >) G ™* h (t', t). 



(21) 



(22) 



evaluated at the saddle-point (here and below the 
summation over repeated indices is assumed). 
The function I((f>) is the mean-field action 



C ij 

-i^TrlnzG- 1 ^) 



Hi 



The expression for M^> is 



M 



fc(i) 



(t) = M* (0) (t) 



+ [6 ls 6 nr S(t 2 - t 3 ) + m{t2M)D%{UM)] 

xBl${tzMM)Dfi{tut), (23) 
where the third order vertex is 



i Sx%(t 3 ,t2) 
2 StfAti) 



(24) 



x (G^™ 1 (ti, t 2 )G^ ri (t 2 , t^G^f 1 (£3, t 2 ) 

+G% Tl ( t i,h)G r s 2 i ni (t 3 , t 2 )G r lf 1 {t 2 ,ti)). 

Let us note that the functions G, \ an d B are di- 
agonal on site indices. The Green's functions G in 
these expressions are calculated self-consistently 
by means of equations l|18|l and (|19J) . 



6 



0.5 



0.25 



-0.25 





A=0.0 




A=0.1 




A=0.2 











5 10 15 20 



0.25 



■0.25 







A=0.0 

A=0.1 

A=0,2 



5 10 15 20 




(t-ti)»c 



Figure 4. Time-dependence of one-spin magneti- 
zation after the 7r/2-rotation of the magnetic field 
Hi = (0,0, H) -> H f = (0,H,0) at different val- 



3. ONE SPIN IN A TIME-DEPENDENT 
MAGNETIC FIELD 

As a first example we consider the problem of 
the behavior of the magnetization of one spin in 
the presence of the homogeneous time dependent 
magnetic field H(t) = (0, 0, H z (t)). We will con- 
sider the effect of performing spin rotations by it 
or it/2. Let us consider first the case of magnetic 
field flipping. It is supposed that the system is in 
equilibrium at initial time t = U at temperature 
Ti and at t = tj + 5 the magnetic field suddenly 
changes its direction: H — > —H. 

Also at this time the coupling to the phonons 
at temperature Tf is switched on, otherwise the 
relaxation to equilibrium will not take place. The 
heat bath is necessary to get magnetization relax- 
ation to equilibrium with the magnetic field —H. 
Otherwise, the magnetization will precess around 
the old field H, as it follows from the equation 
for magnetization in the case of no spin-phonon 
coupling: 

i—M = —M x H. 
dt 

Let us get the equation for 4>f^ (t) when there 
is spin-phonon coupling. Since the magnetic field 
is homogeneous in space, the field 4>[°\t) and the 
magnetization M^{t) do not depend on site i, 
and the mean-field equation (|18|) can be written 
as 

$ 0) (t) =Hi(t)+i [ dt'A(t : t r )M^(t'),(25) 
This last equation is equivalent to 
$f ) {t)=H i {t)- 2 J dt'^a(t-t')MP(t'), (26) 
where 

POO 

%a(t - t') = -A dLuuj s e-" /uJ " sin(w(< - t')) 
Jo 

This is a simple renormalization of the mag- 
netic field induced by the heath bath and does not 
depend on the heat bath temperature Tf . There- 
fore, a higher correction beyond the saddle-point 
level must be considered. 

We shall study the behavior of the magnetiza- 
tion M« of one spin calculated by expressions 
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(|23|l and (|24p. For one spin there are no site in- 
dices in these expressions and the formula for the 
magnetization can be simplified. Let us limit our 
consideration to the linear term in A, which as- 
sumes that the heat bath coupling is weak. Then, 
the expression for the z-component of the mag- 
netization has the following form: 

M z{1 \t) ~ M (0) (i) 



Ci+C 2 



<ftidtaA nn (ti - t 2 )e nrir2 e nnin2 e kklk2 

x [G k2ni (t, h)G n2ri (hM)G r2kl (t 2 , t) 
+G fe2ri (t, t 2 )G r2n > {t 2 M)G n2kl (h, t)]. (27) 

Since we are interested in the correction to the 
magnetization linear in A, the Green's function 
G(ti,t 2 ) in equation f2"7)l can be calculated at 
$°{t) = -H{t). When the magnetic field flip- 
ping is applied, the expressions for free Green's 
functions (|12|1 and l|13l) remain correct with the 
only substitution H z = —H in the complex ex- 
ponents. The Fermi-factors in these expressions 
do not change with the magnetic field flipping, 
since they are defined by the initial conditions, 
but the time exponentials, which define the dy- 
namics of the spins, are changed accordingly. The 
overall expression does not get modified because 
both fields are along the same direction. 

Assuming, for simplicity, that only the x- 
component of the spin is coupled to the heat bath 
(A™" ~ 5 nx ), substituting JTTJ, (O and in 
(|27|l and integrating over t\ and t 2 leads to the 
following expression for the z-component of the 
magnetization. 



M z{1 \t) = M^-A / dujuj s e 



x!l ^M<.w M fe) r (M+HV 

+ 4 + M <°> cotta r ( ^, )2 },(2 8 ) 

where = |tanh(^^) is the initial magne- 

tization. Below we shall use an approximation 



- to) in the numerical calculations, 



for convenience. Obviously, expression H28J1 is cor- 
rect for small times (t — ti), when the perturba- 
tion theory is valid, which is the regime where the 



magnetic field flipping is important. The time de- 
pendence of the magnetization for different types 
of the phonon heat bath is presented in Fig. 2. 
As it follows from this figure the phonon density 
of states plays an important role in the spin dy- 
namics. The behavior of the magnetization after 
lowering the temperature for the cases of mag- 
netic field flipping and constant magnetic field is 
presented in Fig. 3 for comparison. The magne- 
tization decreases as (t — ti) 2 at short times in 
both cases and is a linear function of t — ti, when 
(t — ti) ;§> 1/\H\ (see iffify and asymptotic expres- 
sion H29|l for long times). 

Let us now consider a qualitative scheme of how 
to obtain from expression l|28fl the relaxation of 
the magnetization to the equilibrium value at T = 
Tf. In the case when At = (t — ti) ^> l/ui c , 1/H Z 
the following relation can be used 



sm 2 ((u;±H z ){t-U)/2) tf 



(t-t,)8(uj±H z ).(29) 



(uj±H z ) 2 4' 
Then, equation (|28() can be transformed to 

M z ^(t l + At)=M z( -°\u) 
-XAt(M z ^(U)-Mf 0) ), (30) 
where 



A = 



ttA\H z 



\M 



j/u 



(31) 



/ 



and the equilibrium value of the magnetization at 
T = T f is 



Mf^tanh^) 



(32) 



In the case when At is not too large, but much 
bigger than l/u) c & n d and taking ti as arbi- 

trary, equation (|30|) can be transformed into the 
differential equation 

j t M z{1 \t) = -\{M z{1 \t) - M z{1) ) (33) 
with the initial boundary condition 
M z ^{t t ) = M z ^ = l -t & nh{^). (34) 
The solution of ([Sljl. is 

M z( U(t) = M z{1) + (M ( ^ Z -M z(1) )e- X{t - U) .(35) 
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As it follows from equation (|35|l . the magneti- 
zation approaches its equilibrium value exponen- 
tially, and the relaxation time 1/ A is defined by 
the final temperature, magnetic field, and by the 
heat bath parameters. 

As another example let us now consider the 
time-dependence of the magnetization under the 
effect of a 7r/2-rotation of the magnetic field: 
H, = (0,0, H) -> H f = (0,H,0). The depen- 
dence calculated by means of expression (I27|l in 
the case of the lowest order spin x-component- 
phonon coupling is presented in Fig. 4. The mag- 
netization is precessing and its direction is chang- 
ing towards the new direction of the magnetic 
field in the case of non-zero spin-phonon coupling. 

4. THE CASE OF MANY SPINS 

In the case of many spins the time relaxation 
dynamics is complicated by the spin-spin interac- 
tion. Let us simplify the problem by taking into 
account the spin-spin interaction at the mean- 
field level. The spin Hamiltonian in the case of 
the magnetic field in z-direction can be written 
as 

H spin (t) = -2dJM^S z (t) - H z (t)S z {t), 

where d is dimensionality of the system and the 
initial magnetization is defined by the equa- 
tion 



2 ^ 2 



(36) 



Therefore, the problem is reduced to one spin 
in the effective field H = H z + 2dJM ( -°\ The 
formula (|28|l with the substitution H* — > H z + 
2dJM^ describes the time-dependence of the 
magnetization of the Heisenberg model after the 
magnetic field flipping. The time-dependence of 
the magnetization for the spin system at differ- 
ent values of dJ is presented in Fig. 5. As it fol- 
lows from this Figure, the time-dependence of the 
magnetization strongly depends on the relation 
\H\/dJ keeping the other parameters fixed. Let 
us analyse the expression at quite long times, 
when the relation (|29|l can be applied. Either 
the first or the second terms in the curled brack- 
ets contribute to the integral in H28|) . when the 




Figure 5. Time-dependence of the magnetization 
for spin system after magnetic field flipping at 
different values of dJ (in units of oj c ) and A = 
0.025, s = Q.b,H/uj c = 0.1,T;/w c = 0.1,T f /w c = 
0.02. 



sign of H = —\H\ + 2dJM^ is negative or pos- 
itive, respectively. This leads to different, nega- 
tive or positive, corrections to the magnetization 
due to the magnetic field flipping. Therefore, the 
equation H = defines the critical magnetic field 
value Her-, below which the magnetic field flipping 
does not lead to a change of the magnetization 
sign. Equation H = together with Ij36|l gives 
the expression, which connects H cr with J and 

-di = tanh ^)- 

The initial temperature dependence of H cr is pre- 
sented in Fig. 6. As it follows from this Figure, 
at low temperatures, the critical value of mag- 
netic field is dJ, i.e. it is of order of the ferro- 
magnetic coupling energy, which should be over- 
come to change the sign of the magnetization of 
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Figure 6. Temperature-dependence of H cr . 
Mean-field critical temperature of ferromagnetic 
transition is T c = O.bdJ. 



the system. The time-dependence of the mag- 
netization for the spin system at different values 
of H is presented in Fig. 7. As it follows from 
this Figure, the magnetization does not flip when 
the magnetic field is lower than the critical value 
H cr ~ 0.9575<iJ — 0.1915a-> c at given temperature 
Ti = 0.1o; c (see Fig. 6). In this case the magnetiza- 
tion oscillates with a reasonably large amplitude, 
similarly to the case of the </> 4 -theory 32 33 : and 
does not relax to the overturned magnetization. 
In the 4 -theory a vector self-interacting theory 
is considered. In our case the interaction part is 
due to the spin interactions. 

Finally, in Fig. 8 we present the normalized 
magnetization deviation (M^ - M (£))/[(* - 
ti)u> c ]. The rate of change of the magnetization 
corresponds to the Fermi golden rule for the scat- 
tering cross section. After an initial rapid increase 
the rate converges to a constant valid for interme- 
diate times, as well known. For very large times 



Figure 7. Time-dependence of magnetization for 
spin system at different values of H (in units of 
u) c ) and A = 0.025, dJ/u c = 0.2, s = 0.5, Ti/u c = 
0.1, Tf/u> c — 0.02. Critical value of the magnetic 
field in this case is H cr ~ 0.1915w c . 



the perturbation theory is no longer valid. In 
this regime we have to consider higher order cor- 
rections which complicates the solution. This will 
be considered in a future publication. 

5. CONCLUSION 

The problem of the time-dependence of the 
magnetization of a spin system in a time- 
dependent magnetic field is interesting both from 
the point of view of its practical applications, and 
from the theoretical point of view. In this paper 
we have studied the behavior of the Heisenberg 
ferromagnet coupled to a heat bath in a time de- 
pendent external magnetic field. In particular, we 
considered the case when its direction is suddenly 
changed. As expected the behavior of the system 
strongly depends on the phonon bath properties. 

It has also been shown that there is a critical 
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Figure 8. Time-dependence of the rate of change 
of the magnetization for the spin system after 
magnetic field flipping at different values of s and 
A = 0.025, dJ/uj c = 0.05,H/uj c = 0.1,T,/cj c = 
0.1,T f /u c = 0.02. 



value of the magnetic field below which the mag- 
netization does not relax to the equilibrium value 
after the magnetic field flipping. This situation is 
analogous to the </) 4 -theory case, where also exists 
a critical value of the external source 32 33 . In 
this work we considered a microscopic description 
for the vector spin system directly instead of using 
an effective bosonic theory from the start. Repre- 
senting the spin operators by Majorana fermions 
we constructed the path integral representation 
for the spin system. Also, we used a physical way 
to introduce the relaxation to equilibrium cou- 
pling the system to a heat bath. For simplicity 
we used a phonon heat bath. The spin-phonon 
coupling can be integrated exactly introducing 
an effective time-dependent interaction between 
the spins. Afterwards we bosonized the model 
via a Hubbard-Stratonovich transformation lead- 



ing to an effective bosonic theory (after integrat- 
ing out the Majorana fermions) which has some 
resemblances with the 4> -theory but with signif- 
icant differences like time-dependent coefficients 
and the introduction of cubic terms (like <fi 3 ) due 
to the presence of a three field vertex which origi- 
nates in the three-spin correlation function which 
does not vanish in general due to the nature of 
the spin commutation relations. 

At the same time there are several questions 
which remain open like the interpolation between 
the short time and the long time behaviors. Also, 
in the problem of quantum spinodal decomposi- 
tion, the possibility of the formation of a mag- 
netic bubble of the stable phase immersed in a 
background of the unstable phase and the de- 
termination of its radius time-dependence is an 
interesting open problem. These and some other 
questions are planned to be studied in the nearest 
future where a detailed comparison of the Heisen- 
berg model with the -theory will be carried out. 
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